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ABSTRACT 

We present time-series B, V photometry of NGC 6229, obtained with the purpose of 
providing the first extensive CCD variability study of this cluster. As a result, we were 
able to obtain periods for all NGC 6229 variables, with the exception of five stars from 
the Borissova et al. (1997) candidate list located very close to the cluster center. Two 
variable stars suspected by Carney, Fullton, &; Trammell (1991) are first-overtone RR 
Lyrae (RRc) stars, whereas seven of the 12 candidates from Borissova et al. (1997) are 
confirmed variables — three of them being fundamental RR Lyrae (RRab) pulsators, 
two first-overtone pulsators, one eclipsing binary (most likely an Algol system) and 
one bright star whose variability status could not be satisfactorily determined. A new 
image subtraction method (ISM) suggested by Alard (1999) has been employed which, 
together with the Welch & Stetson (1993) technique, has allowed us to discover twelve 
new RR Lyrae variables in the cluster, for which ephemerides are provided. Ten of 
these are RRab's, whereas the other two are RRc's. As originally suggested by Mayer 
(1961), NGC 6229 is clearly an Oosterhoff type I globular cluster. We also confirm 
that V8 is a Population II Ccphcid of the W Virginis type, and suspect that this is 
the case for V22 as well. 

The physical properties of the NGC 6229 RR Lyrae population are compared 
against those of M3 (NGC 5272) using several different methods, including a standard 
period-shift analysis. Possible differences between these two clusters are discussed. 

Key words: stars: fundamental parameters - stars: horizontal branch - stars: Pop- 
ulation II - stars: variables: other - Galaxy: globular clusters: individual: NGC 6229 



1 INTRODUCTION 

In previous papers of this series, we have carried out exten- 
sive analyses of the colour-magnitude diagram (CMD) of the 
outer-halo globular cluster (GC) NGC 6229 (C1645+476). 
Borissova et al. (1997, hereafter BCSS97) presented the 
first CMD for the crowded core of the cluster, finding the 
horizontal-branch (HB) morphology of NGC 6229 to be bi- 
modal, containing a prominent blue tail with at least one 
'gap'. Catelan et al. (1998) carried out a detailed theoreti- 
cal and statistical analysis of HB bimodality and 'gaps' (in 
NGC 6229 and other GCs). Borissova et al. (1999) obtained 
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the first deep CMD of NGC 6229 clearly identifying the 
main-sequence turnoff of the cluster, and found its metallic- 
ity and age to be very similar to those of M5 (NGC 5904); 
the possible presence of 'extreme' HB stars was also high- 
lighted. 

The purpose of the present paper, the fourth in our 
NGC 6229 series, is to provide the first detailed CCD inves- 
tigation of the variable star population in this very distant 
GC. 

Using unpublished data obtained by H. Shapley, Davis 
(1917) reported on the discovery of the first variable star in 
NGC 6229 — the variable currently known as V8. She called 
attention to the brightness of the variable in comparison 
with other non-variable stars in the cluster, claiming it to be 
the second brightest NGC 6229 star when at its maximum. 
Indeed, Baade (1945) would later rediscover this star and as- 
sign it a 'long-period Cepheid' status. Baade's was the first 
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extensive variability study to be published for NGC 6229. It 
was based on 46 photographic plates obtained in the years 
1932-1935 at the Mount Wilson Observatory. By intercom- 
parison of 11 pairs of plates with a blink-comparator, he 
discovered 20 new variable stars, all 'of the cluster type' 
(i.e., RR Lyrae stars). Sawyer (1953) later extended the list 
of known variable stars in NGC 6229 to 22 with the report 
of a (bright) variable star (V22) very close to the cluster 
core. 

None of these papers reported periods for the variables. 
Mayer (1961) used Baade's (1945) data to provide periods 
for 10 of the RR Lyrae variables. Mannino (1960) published 
extensive photometry for 16 among the 22 known NGC 6229 
variables — and, using both his 92 plates (from Asiago, ob- 
tained over the years 1956-1957-1958) and Baade's 46 plates, 
provided periods for 12 of these stars, assuming that the pe- 
riods stayed constant between the epochs of his and Baade's 
observations. We note that, as a rule, in Sawyer Hogg's 
(1973) catalogue Mannino's periods were adopted, those 
from Mayer being used only for variables V6, V7, and V17, 
for which Mannino did not give periods: for V6 and V7, 
Mannino did not obtain any new data, whereas for V17, 
and also V3, V4, and V21, he was simply unable to derive 
representative periods. Therefore, Sawyer Hogg's catalogue, 
and also its recent update (Clement 1997), contain periods 
for all 22 variables except V3, V4, Vll, V12, V18, V21, and 
V22. 

More recently, Carney, Fullton, & Trammell (1991, 
hereafter CFT91), in their CCD analysis of the CMD of 
NGC 6229, provided single-frame magnitudes and colours 
for most of Baade's (1945) variables. In addition to this, 
they identified two additional stars as likely variables (see 
also Searle & Zinn 1978), besides four other stars for which 
there was some evidence for photometric variability. Finally, 
BCSS97, using the Kadla & Gerashchenko (1982) method 
and a wavelet transform technique to allow CCD photome- 
try in the very dense core of this GC, provided a list con- 
taining an additional 12 RR Lyrae variable candidates. 

In the present paper, we will provide the first CCD light 
curves for all previously known or suspected/confirmed vari- 
ables in NGC 6229.0 Besides, using a new and very power- 
ful 'image subtraction method' (ISM) developed by Alard 
(1999) — the ISIS package — we will also provide light curves 
for twelve new variables discovered in the course of the 
present investigation. 

We begin in the next section by describing our obser- 
vations and data reduction techniques. In Section 3, the 
light curves and the variables' main parameters (periods, 
colours, magnitudes, amplitudes), as obtained using both 
daophot and ISIS, are provided. In Section 4, Blazhko stars 
are identified. In Section 5, the Bailey (period-amplitude) 
diagram is discussed. In Section 6, Fourier decomposition of 
the light curves is used to derive the stellar physical param- 
eters, following methods developed by Simon and Clement 
(RRc stars) and by Kovacs and Jurcsik (RRab stars). In 
Section 7, the nature of a few particularly noteworthy clus- 



t For V3, which is outside our observed field, we shall determine 
the period and variability class on the basis of the data published 
by Mannino (1960). 



ter variables is discussed. Finally, our concluding remarks 
are presented in Section 8. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations 

Time-series photometry was obtained during the interval 
1997 June 28 - 1999 September 26 with the 2m Ritchey- 
Chretien telescope of the Bulgarian National Astronomi- 
cal Observatory 'Rozhen' with a Photometries 1024 x 1024 
CCD camera. The scale at the Cassegrain-focus CCD was 
0.33" pixel -1 and the observing area was 5.6' x 5.6', centered 
on the cluster centre. 

Approximately 200 B and V frames were obtained, with 
exposure times ranging from 120 to 300 sec. The seeing dur- 
ing these observations was between ~ 1" and 1.6". 

2.2 Data reduction: daophot 

All frames were reduced using daophot (Stetson 1993) and 
transformed to the standard system independently. In the 
cases when the seeing was not stable the data were dis- 
carded. The data reduction and absolute photometric cal- 
ibration of the frames are the same as those obtained and 
used for the construction of the deep CMD of NGC 6229 
(Borissova et al. 1999). The internal error of each mea- 
surement for the stars in NGC 6229 depends on the see- 
ing and crowding conditions and is on average w 0.02 mag 
for V < 21 mag and B < 21 mag. For more details 
on the error and completeness analysis, see BCSS97 and 
Borissova et al. (1999). A comparison with the photome- 
try of CFT91 for the non-variable stars in the magnitude 
interval B, V = 15 — 20 mag shows that the two datasets 
are in good agreement, although in some frames we can- 
not exclude systematic zero point differences, though not 
larger than 0.02-0.03 mag. When possible, datapoints from 
the photometry lists of CFT91 and BCSS97 were included 
in the present investigation. 

In order to search for new variables in the studied 
field of NGC 6229, we applied the Welch & Stetson (1993) 
method. Since we have observations made nearly simulta- 
neously in two passbands (B and V) and in several differ- 
ent epochs separated at least by a month, we can calculate 
the Welch & Stetson variability index '/' for all stars. As 
a result, two previously unknown variable stars — V37 and 
V38 — were detected. They are both fundamental RR Lyrae 
pulsators (see below). 

2.3 Data reduction: ISIS 

Alard (1999) and Alard & Lupton (1998) presented a 
new ISM ('isis') aimed at obtaining better light curves in 
crowded fields such as GC cores. The capabilities of the 
method have recently been demonstrated by Olech et al. 
(1999) in a study of variable stars in the core of M5. 

In order to improve the quality of the light curves of 
the known variables in NGC 6229 and to search for new 
variables, all our CCD frames were re-reduced using Alard's 
(1999) code. Ours is one of the very first applications of the 
method in a study of stellar variability in GCs. 
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Twenty-nine of our RR Lyrae stars have magnitudes 
measured using both the standard daophot package (which 
places them in the standard Johnson system) and the non- 
standard ISIS package (which is only capable of providing 
relative fluxes, measured in ADU) . We have used the B and 
V magnitudes for such stars in order to try to obtain trans- 
formation equations linking standard magnitudes and ISIS 
fluxes. As far as we are aware, such a discussion has not 
been presented in the literature so far, the ISIS results usu- 
ally being reported in ADU only (see, e.g., Olech et al. 1999 
in the case of M5) . 

Unfortunately, for most of our stars, the transforma- 
tions appear to be affected by some non-linearity. More 
specifically, the brightest point in the light curves, as ob- 
tained in the standard system, usually appears fainter than 
that calculated from the ISM fluxes. If, as a first approxi- 
mation, we use a linear transformation of the form miSM = 
a + P x mf^ 3 (where m represents ISM magnitudes in B or 
V), the f3 coefficient turns out to be 0.00007 and 0.00014 for 
the V and B filters, respectively, and appears to be fairly 
uniform for our RR Lyrae sample. The zero-point coeffi- 
cients a, however, do seem to vary slightly from star to star, 
so that we cannot convert the ISM fluxes into their Johnson- 
equivalent magnitudes in an entirely consistent way. In the 
present paper, for the RR Lyrae stars which were measured 
only with ISIS, approximate ISM-based B and V magnitudes 
were obtained using the /3 value stated above; the adopted 
a was taken from the standard magnitude measured in the 
reference frames. 

These problems notwithstanding, our analysis clearly 
indicates that ISIS is an extremely powerful tool to de- 
tect new variables and determine their periods, particu- 
larly in very crowded fields. Indeed, applying the method to 
NGC 6229, all known variables in our observed field (which 
excludes only V3) were rediscovered, and light curves of as 
many as ten previously unknown variables were obtained. 
We label these new discoveries V39-V48 in what follows. 



3 PERIODS AND LIGHT CURVES 

Our 'target list' contains 50 variable and possible vari- 
able stars based on Baade (1945), Sawyer (1953), CFT91, 
BCSS97 and the present work (see Sect. 2). We have re- 
analysed the periods of all Baade stars and determined the 
periods of suspected ones using a least-squares periodogram 
analysis by means of the phase dispersion minimization 
(PDM) task available in iraf and a period-finding program 
based on Lafler & Kinman's (1965) 'theta' statistic. We ob- 
tained new light curves and re-derived the periods for 20 
variable stars listed in Sawyer Hogg's (1973) catalogue, and 
also the first light curves for RR Lyrae candidates No. 155 
and No. 88 from CFT91 and 3, 7, 8, 9, 10, 11 and 12 from 
BCSS97. We have not detected any significant variability 
(at a level ^ 3a) for the stars No. 4, No. 105 and No. 134 
suspected by CFT91. A finding chart for the NGC 6229 vari- 
ables is provided in Fig. |l] 

The B and V RR Lyrae light curves obtained with 
daophot are displayed in Fig. and Fig. [], respectively. 

Using the same statistical approach — the PDM task and 
Lafler-Kinman's 'theta' statistic — we analyzed relative ISM 
fluxes of all variables. In general, we recovered the same 



periods obtained with daophot, but in some cases we were 
able to improve the derived periods significantly. The light 
curves in ISM relative fluxes, as obtained from the B and 
V data for the RR Lyrae stars using ISIS, are displayed in 
Fig. ^ and Fig. [B|, respectively. 

The quality of the light curves, particularly for the vari- 
able stars in the inner regions of NGC 6229, is noticeably 
better in the ISIS than in the daophot case (compare, e.g., 
V36 in Fig. ^ and Fig. ^). However, even with ISIS we could 
not obtain representative light curves for candidates num- 
bered 1, 2, 4, 5 and 6 by BCSS97. All these stars are less 
than 0.3 arcmin from the cluster centre (see Fig. |l|) and can- 
not be resolved on our moderate quality CCD material for 
this very distant and fairly dense GC. However, with the 
exception of BCSS97-1, a few of our best ISM-subtracted 
B images do show signs of photometric variability for these 
stars, which thus remain likely variable candidates. Follow- 
up studies are clearly needed to verify these stars' variability 
status. 

In Figs. |2|-|5|, the variable stars are plotted in order of 
increasing period to illustrate the changes in the light curves 
with period. The difference in light curve shape for stars 
with periods greater than 0?40 compared with those with 
shorter period is obvious. Nine variables (V15, V16, V17, 
V23, V24, V32, V36, V44, V47) can be classified as RRc 
stars (i.e., first overtone pulsators). The light curves of the 
remaining 29 variables are more asymmetric and have larger 
amplitudes, and these are undoubtedly RRab Lyrae stars. 

The adopted ephemerides are listed in Tables |IJ |2| ^ 
and[| 

Table |l| contains the list of variable stars discovered by 
Davis (1917), Baade (1945), and Sawyer (1953). Column (1) 
lists the star number from Sawyer Hogg's (1973) catalogue, 
and column (2) gives the variability type of the star. Stars 
identified as Blazhko variables are noted. In columns (3) 
and (4), we give the X and Y coordinates (in arcsec) in the 
Sawyer Hogg (1973) system. Column (5) provides the period 
in days. Columns (6) and (7) give the intensity-mean mag- 
nitudes in B and V, respectively; those values were obtained 
by directly averaging over the pulsation cycle in intensity. 
Likewise, columns (8) and (9) provide the intensity-mean 
colours (B — V) (obtained by directly averaging over the 
colour curves in intensity units) and (B) — {V) [simply ob- 
tained by subtracting the values in column (7) from those 
in column (6)]. Column (10) gives the magnitude-weighted 
mean color, (B — U) ma g, of the stars. Finally, columns (11) 
and (12) provide the estimated amplitudes in B and V, re- 
spectively. 

In order to derive the mean magnitudes and colours, 
the B and V light curves were fitted by smooth curves (e.g., 
Carretta et al. 1998). The smoothing was performed on the 
light curves in intensity units using polynomials of various 
orders. For each star, the ISM light curves were also taken 
in consideration, but only in the sense that all points in the 
B and V light curves which clearly deviate from their corre- 
sponding ones in the ISM light curves were simply rejected. 

We employed Layden's (1998) 'light-curve templates' 
technique to assist us in deriving the RRab amplitudes. 
Lay den only provided templates in the V band. In order 
to derive the B amplitudes, we utilized his V templates as 
well, given the similarity in shape of the normalized light 
curves. We tested the accuracy of the procedure by apply- 
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Figure 1. Finding chart for all known variables in NGC 6229 in a 5.6' X 5.6' area, centered on the cluster centre. North is up and east 
to the left. V3 is outside our field of view. 



ing the V templates to a sample of RRab stars with very 
well-defined B light curves from the literature, encompassing 
a wide variety of light curve shapes. The stars included are 
V9 in 47 Tuc; V3, V7, V10 and V13 in M2 (NGC 7089); V8 
and V28 in M5; V4, V5 and V8 in M92 (NGC 6341); and 
the field star SW Draconis. The data were obtained from 
the following sources. 47 Tuc: Carney, Storm, & Williams 
(1993); M2: Lee & Carney (1999); M5: Storm, Carney, & 
Beck (1991); M92: Carney et al. (1992); SW Dra: Jones 
et al. (1987). As a rule, we found that the B amplitudes 
of these stars were successfully reproduced using the nor- 
malized V templates from Layden. The mean difference be- 
tween the values tabulated in the quoted papers and those 
derived using the V template approach is essentially negli- 
gible, (AAb) = +0.002 ± 0.017 mag (standard deviation). 
However, we recommend that the V templates are not used 
to infer mean magnitudes and colours, since there is a notice- 
able difference in shape between normalized B and V light 
curves at intermediate phases, <f> ~ 0.5. In Figs. 2 and 3, 



we overplot the best-fitting light curves obtained from Lay- 
den's templates on data for the RRab stars with 'regular' 
light curves (see below). 

The entries in Tables 0, ^| and ^ have essentially the 
same meaning as in Table hTtliough the following should be 
noted. In Tables ^ and [| new 'V numbers, starting from the 
last entry in Sawyer Hogg's (1973) catalogue, are assigned in 
order of discovery. For ease of reference, cross-identification 
with the numbers given in CFT91 and BCSS97 is also pro- 
vided. Not listed are the BCSS97 candidates 1, 2, 4, 5, and 
6; however, for future reference, we assign to these likely 
variable stars the numbers V25, V26, V28, V29, and V30, 
respectively. 

Specifically in regard to the new variable stars (V39- 
V48) that we have identified using ISIS (Table^), estimates of 
amplitudes inferred from the minimum and maximum values 
of the relative fluxes are given (see Sect. 2.3 for a discussion). 
For these stars, 'mean' B and V magnitudes and (B) — 
(V) colours are given for reference purposes only, and are 
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Figure 2. B light curves for NGC 6229 variables, as obtained using DAOPHOT. The 'best-fitting' fiducial light curves, as obtained from 
Laydcn (1998), are overplotted on the data for the 'regular' variables (see text). 



based solely on our best B, V image pair. All these uncertain 
values are marked with a colon in Table ^| and must not be 
used except with the purpose of assessing the variability 
type. In any case, it is interesting to note that the RRab 
stars discovered using ISIS tend to have smaller amplitudes 
and longer periods than the bulk of the NGC 6229 RRab 
variables. 

To check the above derived ephemerides, the intensity- 
mean V magnitudes as a function of (B — VOmag are over- 
plotted on the CMD of NGC 6229 in Fig. |. As pointed 
out by Preston (1961) and Sandage (1990), (B - V% ag is 
the most appropriate colour to characterise the temperature 
of an 'equivalent static star'. The NGC 6229 ridgelines are 
from Borissova et al. (1999); for the red HB, it represents 
the lower envelope of the distribution. 

Filled squares denote the RRc stars, and open circles 
correspond to RRab stars. The filled circles represent the 
known Population II Cepheid in the cluster (V8), as well 



as V27 and V31 — whose variability status we will discuss in 
more detail below. 

The mean B and V magnitudes for all RR Lyrae stars 
(except V32, which may be a field variable, and in any case 
is an outlier; see below) for which we have reliable DAOPHOT 
light curves are 18.38 ± 0.10 mag and 18.03 ± 0.08 mag, 
respectively — where the standard deviation of the mean is 
given. The latter value is in good agreement with the HB 
level (Vhb = 18.10 ± 0.05 mag) obtained by Borissova et al. 
(1999; see their Fig. 6) as an estimate of the lower envelope of 
the non-variable red HB stars' distribution in the NGC 6229 
CMD. The mean colour of the RRc stars (except V32 and 
the ISM stars for which we do not have reliable DAOPHOT 
photometry) is {(B — V) mag )c = 0.26 ± 0.03 mag, whereas 
((£> - V) mag )ab = 0.41 ± 0.03 mag for the RRab's. 

Figure ^ suggests that there is some overlap between the 
RRab's and non-variable red HB stars in NGC 6229. How- 
ever, part of this effect appears to be due to the fact that 
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Figure 3. V light curves for NGC 6229 variables, as obtained using DAOPHOT. The 'best-fitting' fiducial light curves, as obtained from 
Laydcn (1998), are overplotted on the data for the 'regular' variables (see text). 



the red HB 'ridgeline' in Fig. 6 of Borissova et al. (1999) 
was indeed extended towards bluer colours than would be 
appropriate for a detailed comparison with RR Lyrae colours 
such as the one we are presently carrying out; a star- by-star 
analysis shows that the overlap is indeed relatively small. 
In order to determine the blue and red edges of the in- 
stability strip, we followed a procedure similar to that of 
Walker (1994). Accordingly, the blue edge of the instability 
strip comes from an average between the two bluest vari- 
ables and the two reddest non- variables on the blue HB, 
giving (B — V)o, be — 0.22 mag. The red edge is simi- 
larly derived by averaging over the five reddest variables 
and the five bluest among the non variables on the red 
HB, which gives (B — V)o, re — 0.45 mag. Also, a 'bound- 
ary' between the RRc and RRab variables can be placed 
at (B — V)o ~ 0.32 mag; there is little, if any, overlap be- 
tween RRc's and RRab's in NGC 6229. A reddening value 
E{B - V) = 0.01 mag was assumed (Harris 1996). 



The periods of RRc stars are between d 264 and d 398, 
with a mean value 0. 32. The periods of RRab stars are be- 
tween d 466 and d 640, with a mean value d 55. In addition, 
the fraction of c-type RR Lyrae variables is f c ~ 0.24. All 
these values are in excellent agreement with NGC 6229's 
classification as an Oosterhoff type I (Ool) GC, as first sug- 
gested by Mayer (1961; see also Castellani & Quarta 1987). 



4 'ANOMALOUS' STARS, BLAZHKO 
VARIABLES AND PERIOD CHANGES 

RR Lyrae stars in GCs, particularly those pulsating in the 
fundamental mode (i.e., RRab's), often show non-repeating 
light curves, which is characteristic of the Blazhko effect 
(e.g., Smith 1995). To identify possible Blazhko variables, we 
used the method derived by Jurcsik & Kovacs (1996), which 
allows one to identify variables with 'anomalous' light curves 
(often because of the Blazhko effect) from analysis of their 
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Figure 4. ISM relative fluxes in B, as derived using ISIS, as a function of phase for the NGC 6229 variables. 



Fourier decomposition parameters. In this method, a 'devi- 
ation parameter' D m is introduced which measures the ex- 
tent to which any given Fourier parameter 'predicted' from 
empirical interrelations among the several other Fourier pa- 
rameters (such relations being based on a large calibrating 
sample of RRab stars, mostly in the field) match the actual 
value of the Fourier parameter for the star. Jurcsik & Kovacs 
argue that if the deviation in any of the Fourier parameters 
can be characterised by D m > 3, the star has a 'peculiar' 



light curve. More recently, Walker & Kovacs (2000) have 
favoured a more restrictive criterion, D m > 5. 

Fourier decompositions of the V light curves were car- 
ried out following Simon & Teays (1982). Fourier decom- 
position amplitudes Ao to Aio and phase differences 4>2i to 
061 were explicitly calculated for use with the Kovacs & 
Kanbur (1998) relations. This was done independently us- 
ing both standard V magnitudes (derived from daophot) 
or ISM relative fluxes (derived from ISIS; see Sect. 2). We ex- 
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Figure 5. ISM relative fluxes in V, as derived using ISIS, as a function of phase for the NGC 6229 variables. 



eluded from our analysis stars which have prominent 'gaps' 
in phase at the maximum of the light curve — cf. V10 and 
V39 in Figs. 3 and 4. We note that the method has been 
calibrated using light curves derived for standard V magni- 
tudes only (Jurcsik & Kovacs 1996; Kovacs & Kanbur 1998), 
so that the D m approach is not strictly applicable to iden- 
tify 'anomalous' stars using our ISM light curves. However, 
since the Fourier phases and the Fourier amplitudes (except 
Ao) do not appear to vary dramatically from one case to the 



other — in spite of some small differences — we shall assume 
that the calibrating equations remain valid in the case of 
iSlS-based 'V magnitudes' as well. 

From this analysis we identified 15 stars with D m > 5 
which can accordingly be classified as variables with 'pe- 
culiar' light curves: V4, V5, V6, V7, V18, V20, V33, V34, 
V35, V37, V38, V40, V41, V46 and V48. For some of these 
stars, as previously stated, light curves are available from 
Baade (1945) and Mannino (1960), thus allowing us to di- 
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Table 1. Periods and mean magnitudes for variables from Sawyer Hogg (1973) 
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a Data from Mannino (1960). 
b Blazhko variable (see text). 
c See Sect. 7.2. 

Table 2. Periods and mean magnitudes for variables suspected by CFT91 



Name 


CFT91 


Type 


X 
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Period 


(«> 


(V) 


(B-V) 


(B) - (V) 


(B — V)mag 


A B 


Ay 




name 
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V23 


No. 88 


RRc 


-19.8 


+57.8 


0.397702 


18.23 


17.91 


0.33 


0.32 


0.31 


0.47 


0.42 


V24 
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RRc 


-12.9 


-58.3 


0.3034368 


18.33 


18.07 


0.27 


0.26 


0.27 


0.60 


0.46 



rectly check for the possible presence of amplitude and/or of photometric errors, period determination and standard fit 

period changes. To achieve this goal, Baade's photographic of the light curves is not larger than 0.05 — 0.07 mag; hence 

(ph) magnitudes were transformed to the Johnson B filter the B amplitude changes appear to be significant at > 3a. 

using 15 comparison stars from Baade's Table 6, whereas Interestingly, all amplitudes are found to be increasing. 

Mannino's were placed by Mannino himself on Baade's sys- For classification purposes, in Tables 1-4 we have 

tern. The periods of these variables were re-derived following assumed that all RRab variables with D m > 5 are 

the same approach as for our CCD data (see Sect. 3 above). Blazhko variables (cf. Jurcsik & Kovacs 1996, especially 

Comparison between the periods that we derived from their Sect. 4.2). However, the possibility cannot be excluded 

Baade's (1945) and Mannino's (1960) data and those listed that there are Blazhko stars whose D m values never become 

in Sawyer Hogg's (1973) catalogue reveals differences which larger than ~ 3, especially if the amplitude modulation is 

are generally less than C ?00009. Since the estimated errors not too strong (Jurcsik 2000). In fact, we note that the RRab 

of our period determination are ^ d 00005, we will consider Lyrae V34 in M3 is a Blazhko variable (Carretta et al. 1998; 

these differences within the (combined) errors and thus con- Clement 1997 and references therein) , whereas Kaluzny et 

sider the periods listed by Sawyer Hogg for these 15 stars al. (1998) have recently obtained D m — 1.9 from their light 

as being satisfactory at > 3a (error of period determina- curve for this star. Conversely, it is well known that there 

tion). As previously stated, these periods were determined are stars with D m > 3 which are not Blazhko variables (see 

by Mannino and by Mayer (1961). Table 8 in Jurcsik & Kovacs 1996). 

Light curves for the 'anomalous' variables V4, V5, V6, In general, comparison between the periods listed in 

V7, V18 and V20 are shown in Fig. [?|. Baade's (1945) data Sawyer Hogg's (1973) catalogue and those derived from our 

are overplotted with plus signs, and Mannino's data with CCD data shows that most of them agree to within d 00009. 

lozenges. Except for V7 — no data available from Mannino Possible significant period changes are found for only two 

(1960) — the amplitude changes for all these stars are obvi- variables: V7 (RRab) and V17 (RRc). For V7 (V17), Sawyer 

ous; they are clearly Blazhko variables. In Table [l|, we quan- Hogg lists a period of P = d 506980 (0 d 324880). Using 

tify the change in B amplitudes as derived from Baade's and Baade's (1945) data we obtain P = d 506978 (0 d 324872), 

Mannino's photographic data and our CCD material. The whereas the period obtained from our CCD data is P — 

error range of amplitude values obtained as a combination d 501600 (0 d 325193). No significant improvement is found 
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Figure 6. Detailed plot showing the RRc (filled squares), RRab (open circles), a Pop. II Cepheid (V8) and two other variables of 
uncertain variability status (filled circles) in the colour-magnitude diagram. The fiducial lines of the cluster (from Borissova et al. 1999) 
are also overplottcd. 



Table 3. Periods and mean magnitudes for variables suspected by BCSS97 
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a Field variable? See Sect. 7.1. 
b Blazhko variable (see text). 



for V17 using Mannino's (1960) data, due to its poor 
quality — Mannino himself was unable to find a period for 
this star. The calculated difference A P p h-CCD for V7 and 
V17 are d 00538 and -0?000313, respectively. Unfortu- 
nately, the reality of these period changes cannot be consid- 
ered definite, since we do not have enough data to construct 
O — C residuals and determine the type and rate of period 
changes. It should be of interest to assemble all the datasets 
available for this cluster (cf. Davis 1917; Baade 1945; Sawyer 
1953; Mannino 1960; and the present work) in order to study 



the period change rates in NGC 6229, given the constraints 
that these can provide on stellar evolution, particularly on 
the HB phase (Smith 1995 and references therein). This will 
be the subject of a future paper. 



5 THE BAILEY DIAGRAM REVISITED 

In Fig. [| the Bailey diagram is plotted. Individual stars are 
identified for ease of reference. Panel a) shows the period- B 
amplitude distribution, whereas panel b) corresponds to the 
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Table 4. Periods and mean magnitudes for variables discovered in the present investigation 
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a Blazhko variable (see text). 



Table 5. B amplitude differences of Blazhko variables 

Name A B A B A B AA B AA fl 

(Baade 1945) (Mannino 1960) (this work) (Baade-this work) (Mannino— this work) 

V4 0.53 0.70 1.15 -0.62 -0.45 

V5 0.70 0.79 1.30 -0.60 -0.51 

V6 0.68 - 1.15 -0.47 

V18 0.66 - 0.90 -0.24 

V20 1.14 1.14 1.60 -0.46 -0.46 



period-V amplitude distribution. RRc variables (short peri- 
ods and low amplitudes) stand out clearly from their RRab 
counterparts in this diagram. NGC 6229 variables from Ta- 
bles 1, 2, and 3 are shown as open circles. Only RRab stars 
with D m < 5 are displayed; the only exception is V10, for 
which we have not computed a D m value (cf. Sect. 4). Com- 
parison between our light curves and those based on the data 
from Baade (1945) and Mannino (1960) does not show any 
sign of amplitude variations at the » 0.1 mag level, which 
suggests that V10 is not a Blazhko variable. 

Also plotted are M3 stars (data from Carretta et al. 
1998) and M15 (NGC 7078) variables (data from Bingham et 
al. 1984, as selected by Catelan 1998). As well known, these 
GCs are typical representatives of the Ool and OoII classes, 
respectively. M3 stars are shown as solid circles, whereas 
open triangles indicate the M15 variables. As was done in 
the case of M15 by Catelan, we have removed all M3 stars 
that Carretta et al. explicitly note to be (possible) Blazhko 
variables, as well as M3 variables identified as showing the 
Blazhko effect in Clement's (1997) catalogue. In addition, we 
cross-checked the Carretta et al. and Kaluzny et al. (1998) 
lists to identify possible stars in common between the two 
studies with large D m values. None was found. As already 
mentioned though, we found that V34, which is a Blazhko 
variable according to Carretta et al. and Clement's cata- 
logue, has a very small D m value — which might have led 
one to conclude that it has a 'normal' light curve. V34 was 
also excluded from Fig. |^. 

It is immediately clear that there is a period shift (at 
fixed amplitude) between NGC 6229 and M3 RRab vari- 
ables, in the sense that M3 stars have longer periods than 
their NGC 6229 counterparts. This is more easily noticed 
when we overplot the lines that best fit the M3 data, ob- 



tained following the 'OLS bisector' method as appropriate 
in this case (cf. Isobe et al. 1990): 

log P ab (M3) = -0.097 - 0.127 As, 

log P ab (M3) = -0.110 - 0.151 Av- 

Most NGC 6229 RRab stars with regular light curves fall 
below these best-fitting M3 lines, indicating shorter RRab 
periods at a given amplitude than is the case for M3. Surpris- 
ingly, the size of the resulting period shifts is not the same 
when the B and V amplitudes are considered. In the former 
case, one finds a mean period shift for the 10 RRab stars 
(in the sense NGC 6229 minus M3) of (A log P(A B )) = 
—0.021 ± 0.009 (standard error of the mean); in the latter, 
one has instead (A log P{A V )) = -0.008 + 0.008. Given the 
slopes of the M3 period-amplitude relations, such a differ- 
ence could be accounted for if the V amplitudes were too 
large by ~ 0.09 mag, or if the B amplitudes were too small 
by ~ 0.11 mag (or some combination thereof). We believe 
that better data will be necessary before the cause of this 
difference can be fully explained. We shall assume that the 
mean between the two values, —0.015, is a good indicator of 
the actual period shift between the two clusters. 

Fundamental-mode blue amplitudes and temperatures 
appear to be closely related (Sandage 1981a, 1981b; Jones 
et al. 1992; Catelan 1998; Sandquist 2000). Therefore, pe- 
riod shifts at constant blue amplitude can be interpreted 
as period shifts at (approximately) constant temperature. 
From the period-mean density relation (van Albada & Baker 
1971), it thus follows that the detected period shift must 
be due to luminosity and/or mass differences between M3 
and NGC 6229 variables. Assuming similar masses for the 
variables in the two clusters, this implies that the M3 vari- 
ables might be brighter than their NGC 6229 counterparts 
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Figure 7. B light curves as a function of phase for Blazhko variables. The plus signs and lozenges represent Baade's (1945) and Mannino's 
(1960) photographic magnitudes, respectively, transformed to the Johnson B filter. 



by ~ 0.05 mag. It is worth noting that such luminosity dif- 
ferences at a fixed metallicity, if indeed present, can only be 
caused by variations in 'second-parameter' candidates other 
than age or mass loss on the red giant branch, with possible 
implications for our understanding of the 'HB bimodality' 
phenomenon (e.g., Catelan, Sweigart, & Borissova 1998). On 
the other hand, the period shift could also, at least in prin- 
ciple, be due to a mass difference between the RR Lyrae 
variables in the two GCs. At a fixed luminosity, a mass dif- 
ference AMk 0.03 Mq, in the sense that M3 variables are 
less massive, would be required in this case. 



We recall, from Borissova et al. (1999), that NGC 6229 
appears to have virtually the same metallicity as M5. On 
the other hand, it is quite well established that M5 is more 
metal-rich than M3 by ~ 0.2 — 0.3 dex (see Jurcsik 1995 and 
references therein). Therefore, contrary to Clement & Shel- 
ton (1999), our results suggest that there is no universal, 
metallicity-independent 'Ool line' in the period-amplitude 
diagram. However, we cannot exclude some systematic ef- 
fect due to our relatively small number of RRab stars with 
regular light curves and/or to uncertainties in the derived 
amplitudes. Better data for NGC 6229 are needed before this 
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Figure 8. Bailey diagrams for the RR Lyrae stars in NGC 6229 (open circles). Panel a) shows the period- B amplitude diagram, while 
panel b) shows the period-V amplitude distribution. For comparison purposes, data for the Ool GC M3 (filled circles) and the Ooll GC 
M15 (filled triangles) arc also displayed. Blazhko and/or variables with large D m values (Jurcsik & Kovacs 1996) were avoided (see text). 
Individual RRab stars in NGC 6229 are indicated by their names. The solid line in each plot corresponds to the 'OLS bisector' best fit 
to the M3 RRab data (see text). 



issue can be settled with certainty. We mention, in passing, 
that Kaluzny et al. (2000) have recently claimed that the 
period-V amplitude diagram is indeed a function of metal- 
licity, based on their M5 analysis. 



6 PHYSICAL PROPERTIES OF RR LYRAE 
STARS FROM THEIR FOURIER 
PARAMETERS 

6.1 RRc stars 

Simon & Clement (1993) used light curves of RRc stars ob- 
tained from hydrodynamic pulsation models to derive equa- 
tions to calculate mass, luminosity, temperature and a 'he- 
lium parameter' as a function of the Fourier phase difference 
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Table 6. Physical parameters for RRc variables in NGC 6229 



Name 


A 


Ai 


#31 


Mass 


log (£/£©) 




y 


^(#31 










(M Q ) 




(K) 






V15 


18.14 


0.249 


3.512 


0.512 


1.62 


7526 


0.30 


0.04 


V16 


18.02 


0.178 


3.077 


0.625 


1.72 


7297 


0.26 


0.26 


V17 


17.96 


0.212 


3.120 


0.620 


1.72 


7294 


0.26 


0.07 


V23 


17.93 


0.175 


4.146 


0.532 


1.75 


7178 


0.26 


0.03 


V24 


18.08 


0.206 


3.239 


0.581 


1.68 


7379 


0.28 


0.04 


V32 


18.21 


0.076 


3.513 


0.602 


1.76 


7188 


0.25 


0.17 



Table 7. Physical parameters for RRc variables from ISM light 
curves 



Name 


Ai 


#31 


Mass 


log (L/L B ) 


T e s 


V 


<*(#3l) 








(Me) 




(K) 






V15 


0.184 


3.499 


0.514 


1.62 


7524 


0.30 


0.03 


V16 


0.140 


3.499 


0.562 


1.70 


7337 


0.28 


0.04 


V17 


0.175 


3.564 


0.555 


1.70 


7335 


0.28 


0.02 


V23 


0.192 


3.512 


0.624 


1.79 


7119 


0.25 


0.04 


V24 


0.158 


3.832 


0.500 


1.65 


7435 


0.30 


0.03 


V32 


0.076 


3.320 


0.637 


1.78 


7158 


0.25 


0.08 


V36 


0.152 


3.649 


0.487 


1.60 


7569 


0.31 


0.04 


V44 


0.168 


3.171 


0.644 


1.76 


7199 


0.26 


0.03 


V47 


0.205 


3.611 


0.555 


1.70 


7312 


0.28 


0.03 



</>3i and period. From their equations (2), (3), (6) and (7), we 
computed these quantities using the Fourier decomposition 
parameters obtained from our CCD light curves. They are 
given in Table ^. Note that the 'helium parameter' y does 
not provide a good description of the helium abundance (Si- 
mon & Clement 1993). 

We exclude from our analysis V36 due to its irregular V 
light curve and hence large #31 error. V16 is another variable 
with an irregular V light curve and relatively large a(<f>3i). 

Excluding V32 (see below), the unweighted mean val- 
ues and standard deviations of the mass M, luminosity 
log (L/L0), effective temperature T e g and 'helium param- 
eter' y are: 0.57 ± 0.05 M©, 1.70 ± 0.05, 7335 ± 129 K 
and 0.27 ± 0.02, respectively. For comparison purposes, we 
find, from the latest VandenBerg et al. (2000) a-enhanced 
zero-age HB (ZAHB) models, at T cff = 7300 K, the fol- 
lowing values: for [Fe/H] ~ -1.4 (NGC 6229 metallicity 
in the Zinn & West 1984 scale; cf. Borissova et al. 1999), 
Mzahb ^ 0.62 Mq, log(L/L )zAHB 1.62; for [Fe/H] ~ 
-1.1 (NGC 6229 metallicity in the Carretta & Gratton 
1997 scale; cf. Borissova et al. 1999), Mzahb ^ 0.59 M©, 
log (L/L©)zahb — 1.61. As can be seen, the Simon & 
Clement (1993) method favours the 'long' Pop. II distance 
scale, providing higher HB luminosities than those predicted 
by VandenBerg et al. (2000). 

The same analysis was carried out using the ISM (isis) 
light curves, and the results are presented in Table [?]. Note 
the remarkable decrease in cr(#3i) for V16 and V32 from Ta- 
ble ^ to Table ^| We were also able to successfully analyse 
V36 with the ISM light curves, obtaining a relatively small 
a (#31). The unweighted mean values and standard devia- 
tions of the mass M, luminosity log (L/L©), effective tem- 
perature T e fs and 'helium parameter' y axe: 0.56 ± 0.05 M©, 
1.69 ± 0.06, 7332 ± 157 K and 0.28 ± 0.02, respectively. As 
can be seen, the ISIS values are in good agreement with those 
obtained using the standard V light curves. 

How does this compare with M3? For this GC, Kaluzny 
et al. (1998) have found the following mean values: (M) = 



Table 8. Physical parameters for RRab variables in NGC 6229 



Name 




A-, 
/ll 


<psi 


#41 


[l-e/iij 


Ivly 




^/fr /ul \ 
















(IS.) 




V 1 


17.934 


0.400 


4.949 


1.182 


— 1.54 


0.72 


6438 


0.09 


V2 


18.078 


0.371 


4.800 


1.189 


-1.58 


0.76 


6431 


0.08 


V9 


18.029 


0.273 


5.037 


1.119 


-1.19 


0.85 


6493 


0.15 


Vll 


18.181 


0.288 


4.994 


1.279 


-1.27 


0.81 


6478 


0.09 


V12 


18.043 


0.369 


4.483 


4.591 


-1.45 


0.78 


6444 


0.10 


V13 


18.106 


0.293 


4.858 


1.373 


-1.46 


0.82 


6410 


0.07 


V14 


17.926 


0.360 


4.688 


0.787 


-1.25 


0.88 


6595 


0.09 


V19 


18.067 


0.452 


4.883 


1.005 


-1.04 


0.84 


6654 


0.08 


V21 


18.025 


0.248 


5.392 


1.328 


-0.83 


0.87 


6355 


0.11 



Table 9. Physical parameters for RRab variables in NGC 6229 
from the ISM light curves 



Name 


[Fc/H] 


My 


T c ft 


<r([Fe/H]) 


VI 


-1.65 


0.73 


6370 


0.09 


V2 


-1.61 


0.77 


6423 


0.08 


V5 


-1.53 


0.81 


6462 


0.13 


V6 


-1.41 


0.79 


6308 


0.12 


V7 


-1.65 


0.83 


6380 


0.16 


V9 


-1.46 


0.81 


6376 


0.15 


Vll 


-1.40 


0.83 


6437 


0.09 


V12 


-1.44 


0.79 


6431 


0.10 


V13 


-1.29 


0.83 


6427 


0.07 


V14 


-1.27 


0.89 


6525 


0.09 


V19 


-1.68 


0.82 


6446 


0.08 


V20 


-1.54 


0.81 


6553 


0.19 


V21 


-0.93 


0.87 


6372 


0.11 


V33 


-1.33 


0.85 


6025 


0.09 


V34 


-1.48 


0.81 


6356 


0.13 


V38 


-1.45 


0.84 


6487 


0.18 


V40 


-0.81 


0.84 


6227 


0.15 


V42 


-1.75 


0.81 


6379 


0.17 


V43 


-1.74 


0.74 


6344 


0.13 


V45 


-0.85 


0.85 


6322 


0.12 



0.59 Mq\ (log (L/L©)) = 1.71; (T e «) = 7315 K. This sug- 
gests masses lower in NGC 6229 than in M3 by ps 0.02 - 
0.03 M© (in the opposite sense to that required to account for 
the detected period shift between the two clusters), and lu- 
minosities lower in NGC 6229 than in M3 by 0.03-0.05 mag 
(as needed to account for the period shift between these clus- 
ters). One should note, however, that this method may need 
further testing and calibration, given the empirical results 
recently presented by Kovacs (1998). 

6.2 RRab stars 

Jurcsik & Kovacs (1996), Kovacs & Jurcsik (1996, 1997) 
and Kovacs & Kanbur (1998) obtained empirical formulae 
relating the stellar metallicities, absolute magnitudes and 
temperatures to Fourier decomposition parameters for RRab 
stars with 'regular' light curves. The only model-dependent 
ingredients in their calibrations are the zero point of the 
HB luminosity scale (adopted from Baade-Wesselink stud- 
ies) and the colour-temperature transformations (obtained 
from static model atmospheres). 

The physical parameters of NGC 6229 RRab Lyrae vari- 
ables obtained from this method are given in Table ^. The 
unweighted mean value (and corresponding standard de- 
viation) of [Fe/H] derived from #31 is —1.26 ± 0.24 dex. 
Note that this value is in the Jurcsik (1995) scale. Like- 
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wise, the mean absolute magnitudes and temperatures are 
(M v ) = 0.82 ± 0.05 mag and (T eff ) = 6478 ± 93 K. 
Again, similar values are obtained from the ISM light curves: 
([Fe/H]) = -1.41 ± 0.28 dex, (My) = 0.81 ± 0.04 mag and 
(T cff ) = 6383 ± 116 K. The faint HB level is a reflection of 
the adoption of the Baade-Wesselink luminosity zero point 
in the calibration of this method (see Jurcsik & Kovacs 1999 
for a recent discussion). 

Comparison with the mean absolute magnitude of the 
RRab Lyrae in M3 recently obtained by Kaluzny et al. 
(1998) using the same approach, {Mv) = 0.78 ±0.01 mag, is 
again suggestive of a slightly fainter HB level in NGC 6229 
than in M3, as required in order to account for a period 
shift between the RRab variables in the two clusters. For 
an analysis of possible problems related to this method, the 
reader is refered to Kollath, Buchler, & Feuchtinger (2000) 
and Jurcsik & Kovacs (1999). 



7 COMMENTS ON INDIVIDUAL STARS 

V3: Although data (of low quality) have been provided for 
this star by Baade (1945) and by Mannino (1960), no period 
has been published to date. We note that our CCD investiga- 
tion does not provide help in resolving the situation, since 
V3 is not in our observed field. However, we have investi- 
gated Mannino's data and found that it was possible to fold 
his observations with a period P — 0? 575218. Fig. ^ shows 
the first light curve for this star, which is clearly an RRab 
Lyrae variable. Due to its distance from the cluster centre 
and poor quality of the light curve, V3 was not considered 
when deriving mean cluster properties. 

V12: CFT91 noted that this RRab star, in comparison with 
the other NGC 6229 variables, 'seems to be much too bright 
and too blue'. They suspected that this could be due to 
crowding, hence blending. 

Careful inspection of our images does reveal that V12 
has a very close 'companion'. This is the main reason why 
the ISM light curves for this star (Figs. ^ and are signifi- 
cantly better than their daophot counterparts (Figs. ^| and 
^|): ISIS encountered no difficulty at removing this close com- 
panion. Note that, unlike in CFT91, the mean magnitudes 
and colours for this star (Table 1) do not deviate substan- 
tially from the remainder of the sample. 



7.1 RR Lyrae variables 

V32: This RRc star is the faintest overtone variable in the 
NGC 6229 CMD (Fig. |). V32, at (V) = 18.27 mag, is almost 
0.2 mag fainter than the lower envelope of the red HB of the 
cluster. On the other hand, the Fourier decomposition anal- 
ysis suggests that V32 is one of the two intrinsically brightest 
RRc variables in the NGC 6229 field (cf. Tables 6 and 7). In 
spite of V32's proximity to the cluster centre (Fig. [l]), these 
properties suggest that V32 may be a non-member. Accord- 
ingly, all average values reported in this paper do not take 
V32 into account. 



7.2 Population II Cepheids 

V8: As previously stated, Davis (1917) found V8 to be one of 
the brightest stars in NGC 6229, and Baade (1945) classified 
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Figure 9. B light curve for V3, obtained from Mannino's (1960) 
ph data transformed to the Johnson system using the indicated 
period. 
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Figure 10. Light curves as a function of phase for V8. The plus 
signs and lozenges represent Baade's (1945) and Mannino's (1960) 
photographic magnitudes, respectively, transformed to the John- 
son B filter as described in the text. 



it as a 'long-period Cepheid'. Mannino (1960) obtained a pe- 
riod P — 14"?845093, thus assigning to it W Virginis status. 
The period we estimate from our new CCD data is slightly 
shorter: P = 14?840457. The period, B amplitude and light 
curves (Fig. |l(]) are all consistent with its classification as 
a W Vir variable. In comparison with the typical colours of 
W Vir stars in GCs, V8, with (B - V)o, mag ^ 0.77 mag, 
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Figure 11. B, V light curves for V22, obtained using the ISM 
method (isis) as applied to our shortest-exposure frames. 



seems to be among the reddest: from Harris (1985), one 
might have expected a colour (B — V)o — 0.52 mag given this 
star's period, whereas from Nemec, Nemec, & Lutz (1994) 
a colour (B — V)o — 0.6 mag would seem more appropriate. 
Using the distance modulus from Harris (1996) and given 
that (V) = 15.53 mag, its absolute magnitude is found to 
be (tfv) — —1.85 mag, which places it among the funda- 
mental Pop. II variables in the Nemec et al. classification 
scheme. 

V22: No light curve has ever been published for V22, which 
is located very close to the core of the cluster according to 
the coordinates provided by Sawyer (1953). To the best of 
our knowledge, no follow-up observations were carried out 
for this star since Sawyer reported on its discovery, although 
Sawyer Hogg's (1953) catalogue mentions that this variable 
is 'probably slow' (i.e., has a long period). 

In order to assess the variability type and attempt to 
determine a period for this star, we have used exclusively our 
10 best and shortest (exposure times 60 sec) ISM-subtracted 
frames, since we expected that any star that could be found 
by Sawyer (1953) at the very core of the cluster using her 
photographic material had to be very bright, and hence 
probably only seen by Sawyer on her frames obtained with 
the shortest exposures. This expectation was also reinforced 
by Sawyer's remark that the star's 'magnitude indicates that 
its period may be greater than a day'. 

The ISM method revealed no variable star at the po- 
sition reported by Sawyer Hogg (1973). However, relatively 
close to her coordinates, at X — +3.8" and Y = —10.5", a 



very bright variable star was found. We suggest that this is 
probably a rediscovery of Sawyer's (1953) V22, and that the 
previous coordinates for this star were incorrect. 

Since we also have ~ 10 daophot magnitudes measured 
on our short (60 sec) frames, we attempted to derive a trans- 
formation between ISM relative fluxes and standard B, V 
magnitudes for this star. We found that it was necessary to 
multiply the V22 ISM fluxes by 0.00002 and 0.0000046 in 
order to obtain magnitudes in B and V, respectively. Note 
that these coefficients are different from those obtained for 
the RR Lyrae variables in Section 2.3. One possible reason 
for such a difference is related to the sky subtraction in this 
extremely crowded region of the cluster (see Fig. 1). Using 
B, V magnitudes thus derived, light curves are shown in 
Fig. [l(]. As can be seen, a period P = 15? 8373 gives an 
acceptable light curve for the star, suggesting that it is a 
W Vir variable. 

The derived photometric parameters for this star are 
very uncertain, but they allow us to constrain the variabil- 
ity status of V22. The colour of the star is (B — V)o ~ 
0.9 mag, making it the reddest among all known NGC 6229 
variables. From Harris (1985), one would expect a colour 
(B — V)o — 0.52 mag given this star's period, assuming it 
to be a W Vir variable, whereas from Nemec et al. (1994) a 
colour (B — V)o — 0.6 mag would seem more appropriate. 
Using the distance modulus from Harris (1996), V22's abso- 
lute magnitude is found to be (My) ~ —2.3 mag. This would 
place V22 among the overtone Pop. II Cepheids, following 
the Nemec et al. classification scheme. Note, however, that 
McNamara (1995) has questioned the reality of the subdi- 
vision of Pop. II Cepheids into fundamental and overtone 
pulsators by Nemec et al. (1994). 

In summary, we confirm that V8 and V22 are most likely 
W Vir variables, but more data are needed, especially in the 
case of V22, to verify their colours. 

7.3 Other variable stars 

Two additional stars, V27 and V31, have average magni- 
tudes placing them above the HB level of NGC 6229 (see 
Fig. p|). These bright variables were both discovered by 
BCSS97. 

V27: V27 has a period P = l d 13827. The light curves of V27 
in the relative ISM B, V fluxes are presented in Fig. [l2|. V27 
has period, colour and luminosity appropriate for its classi- 
fication as a Pop. II Cepheid of the BL Herculis type. How- 
ever, the light curve is not well defined, and we regard the 
derived period uncertain, as well as its BL Her classification. 
V31: From the analysis of daophot magnitudes, V31 has a 
period P = 0'?715208. A different period is derived on the ba- 
sis of the ISM fluxes; in this case, we find P = $698893. V31 
appears to lie on an extension of the BL Her period-colour 
relation (Nemec et al. 1994) towards bluer colours, and its 
luminosity would be consistent with this interpretation. V31 
might also be classified as an 'AHB1' star (Diethelm 1983). 

This having been stated, we emphasize that the light 
curve (Fig. does not obviously resemble either those of 
BL Her or of AHB1 stars (cf. Diethelm 1983). In fact, the 
light curve remarkably resembles that of an eclipsing binary. 
V31 might be a contact binary (W UMa star). We have 
compared the properties of V31 against those of W UMa- 
type binaries in GCs from Rucinski (2000). Judging from 
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Figure 12. ISM light curves for V27, obtained using ISIS. Figure 13. ISM light curves for V31, obtained using ISIS. Super- 

imposed on the data is a light curve 'template' for Algol variables 
(from Layden & Sarajedini 2000). 



Rucinski's eq. (1) (a period-colour-luminosity relation), the 
period and colour of the star are inconsistent with its clas- 
sification as a cluster W UMa star: V31's implied absolute 
magnitude, Mv ~ 1.54 mag, corresponds to an apparent V 
magnitude of V m 19.0 mag, adopting a distance modulus 
(m - M)v = 17.46 mag (Harris 1996). This is ~ 1.6 mag 
fainter than V31, which has (V) — 17.4 mag (cf. Table 3). 
Alternatively, the star might be a field W UMa variable. 
However, V31 is very close to the cluster centre, making it 
an unlikely (but not impossible) field variable. In addition, 
the colour and amplitude of the star, compared against the 
results compiled by Rucinski, appear to militate against the 
very classificaion of V31 as a W UMa star. In particular, it is 
clear from Rucinski's Table 2 that the bluest W UMa stars 
(with colours similar to V31's) tend to have much smaller 
amplitudes than is the case with V31. The star also lies 
outside the domain of W UMa variables as shown in Rucin- 
ski's Fig. 6 (period-colour plane). Therefore, we suggest that 
this star is instead a short-period detached or semi-detached 
(most likely an Algol, but possibly a /3 Lyrae) foreground 
eclipsing binary. This classification is certainly consistent 
with the star's light curve, as shown in Fig. [13|, where the 
unequally deep eclipses stand out, and there is an indication 
of a period during which the light may be constant; accord- 
ingly, the Algol 'template' from Layden & Sarajedini (2000) , 
kindly provided by A. Layden, seems to provide a good fit 
to the star's light curve (Fig. |l3j ). 

However, as in the case of V27, we consider the periods 
and light curves of V31 preliminary. Additional data are 
needed to determine its variability status and ephemerides 
with greater confidence. 



8 SUMMARY 

In the present paper, we have presented the first extensive 
CCD investigation of the variable star population in the 
outer-halo GC NGC 6229. With the exception of five vari- 
able stars from the candidate list presented by BCSS97, we 
were able to derive periods and B, V light curves for 36 pre- 
viously known/suspected variables in the field of the clus- 
ter, as well as for 12 newly discovered variables. Most of 
the new discoveries we owe to an application of the new 
ISM method developed by Alard (1999) and Alard & Lup- 
ton (1998), which does appear to be a very powerful tool to 
search for new variables and determine their periods, espe- 
cially in crowded regions such as the cores of GCs. 

The confirmed variable-star population in NGC 6229 
is mostly comprised of RRab stars, several of which present 
the Blazhko effect. Less numerous, but also well represented, 
are the RRc stars. The mean periods of the fundamental 
RR Lyrae pulsators, as well as their relative incidence in 
comparison with first-overtone RR Lyrae variables, clearly 
classify NGC 6229 as an Ool cluster, as originally suggested 
by Mayer (1961). 

We have found that the period-amplitude distribution 
of NGC 6229 RRab variables with 'regular' light curves may 
deviate from that of M3, a very well-known Ool globular. 
The period shift is in the sense that NGC 6229 variables 
have shorter periods than their M3 counterparts at a given 
amplitude. Given that NGC 6229 appears to have the same 
metallicity as M5 (Borissova et al. 1999), which in turn is 
well known to be ~ 0.2 — 0.3 dex more metal-rich than M3, 
this might indicate that the period-amplitude diagram is 
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not solely a function of Oosterhoff type as recently sug- 
gested, but also of metallicity — a conclusion also reached 
by Kaluzny et al. (2000) in the case of M5. However, we 
caution that the size of the period shift was not consistent 
when B or V amplitudes were employed, so that systematic 
effects cannot be ruled out at this point. Better data will be 
necessary before this problem can be conclusively solved. 

We have employed Fourier decompositions of the RRc 
and RRab light curves in an effort to determine their phys- 
ical parameters (Simon & Clement 1993; Kovacs & Jurcsik 
1996, 1997); the results were compared with those for other 
clusters (as published in the literature), and placed in the 
context of our period-shift discussion. 

Finally, we have found that NGC 6229 may contain up 
to two Pop. II Cepheids of the W Vir type, as well as an 
eclipsing binary (most likely a foreground Algol system) and 
an additional bright variable for which we were unable to 
reliably determine the variability type. 



ACKNOWLEDGEMENTS 

The authors are grateful to C. Alard, B.W. Carney, C. 
Clement, J. Jurcsik, L.K. Fullton, A.C. Layden, A. Olech 
and R.T. Rood for very helpful discussions, information, 
and/or kindly supplied data. This research was supported 
in part by the Bulgarian National Science Foundation grant 
with the Bulgarian Ministry of Education and Sciences. Sup- 
port for M.C. was provided by NASA through Hubble Fel- 
lowship grant HF-01105.01-98A awarded by the Space Tele- 
scope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., for NASA 
under contract NAS 5-26555. 



REFERENCES 

Alard C, 1999, A&A, 343, 10 

Alard G, Lupton R.H., 1998, ApJ, 503, 325 

Baade W., 1945, ApJ, 102, 17 

Bingham E.A., Cacciari C, Dickens R.J., Fusi Pecci F., 1984, 

MNRAS, 209, 765 
Borissova J., Catelan M., Ferraro F.R., Spassova N., Buonanno 

R., Iannicola G., Sweigart A.V., 1999, A&A, 343, 813 
Borissova J., Catelan M., Spassova N., Sweigart A.V., 1997, AJ, 

113, 692 (BCSS97) 
Carney B.W., Fullton L.K., Trammell S.R., 1991, AJ, 101, 1699 

(CFT91) 

Carney B.W., Storm J., Trammell S.R., Jones R.V., 1992, PASP, 
104, 44 

Carney B.W., Storm J., Williams C, 1993, PASP, 105, 294 
Carreta E., Cacciari C, Ferraro F.R., Fusi Pecci F., Tessicini C, 

1998, MNRAS, 298, 1005 
Carreta E., Gratton R.G., 1997, A&AS, 121, 95 
Castellani V., Quarta M.L., 1987, A&AS, 71, 1 
Catelan M., 1998, ApJ, 495, L81 

Catelan M., Borissova J., Sweigart A.V., Spassova N., 1998, ApJ, 
494, 265 

Catelan M., Sweigart A.V., Borissova J., 1998, in A Half Century 
of Stellar Pulsation Interpretations, eds. P.A. Bradley & J. A. 
Guzik (San Francisco: ASP), ASP Conf. Ser., 135, 41 

Cle ment G. 1997. AAS Newsl.. 84. 15 



Diethelm R., 1983, A&A, 124, 108 

Harris H.C., 1985, in Cepheids: Theory and observations, IAU 
Colloq. 82, ed. B.F. Madore (Cambridge: Cambridge Univer- 
sity Press), p. 232 

Harris W.E., 1996, AJ, 112, 1487 

Isobe T., Feigelson E.D., Akritas M.G., Babu G.J., 1990, ApJ, 
364, 104 

Jones R.V., Carney B.W., Latham D.W., Kurucz R.L., 1987, ApJ, 
314, 605 

Jones R.V., Carney B.W., Storm J., Latham D.W., 1992, ApJ, 
386, 646 

Jurcsik J., 1995, AcA, 45, 653 

Jurcsik J., 2000, private communication 

Jurcsik J., Kovacs G., 1996, A&A, 312, 111 

Jurcsik J., Kovacs G., 1999, New AR, 43, 463 

Kadla Z., Gerashchenko A., 1982, Izv. Main A. O. of the Academy 
of Sc. USSR (Pulkovo), 199, 86 

Kaluzny J., Hilditch R.W., Clement C, Rucinski S.M., 1998, MN- 
RAS, 296, 347 

Kaluzny J., Olech A., Thompson I., Pych W., Krzeminsi W., 

Schwarzenberg-Czerny A., 2000, A&AS, 143, 215 
Kollath Z., Buchler J.R., Feuchtinger M., 2000, ApJ, 540, 468 
Kovacs G., 1998, MSAIt, 69, 49 
Kovacs G., Jurcsik J., 1996, ApJ, 466, L17 
Kovacs G., Jurcsik J., 1997, A&A, 322, 218 
Kovacs G., Kanbur S., 1998, MNRAS, 295, 834 
Lafler J., Kinman T., 1965, ApJS, 11, 216 
Lee J.-W., Carney B.W., 1999, AJ, 117, 2868 
Layden A.C, 1998, AJ, 115, 193 
Layden A.C, Sarajedini A., 2000, AJ, 119, 1760 
McNamara D.H., 1995, AJ, 109, 2134 

Mannino C, 1960, Pubbl. delFOss. astr. univ. Bologna, Vol. VII, 
N. 13 

Mayer P., 1961, BAICz, 12, 167 

Nemec J.M., Nemec A.F.L., Lutz T.E., 1994, AJ, 108, 222 
Olech A., Wozniak P.R., Alard C, Kaluzny J., Thompson LB., 

1999, MNRAS, 310, 759 
Preston G.W., 1961, ApJ, 133, 29 
Rucinski S.M., 2000, AJ, 120, 319 
Sandage A., 1981a, ApJ, 244, L23 
Sandage A., 1981b, ApJ, 248, 161 
Sandage A., 1990, ApJ, 350, 603 
Sandquist E.L., 2000, MNRAS, 313, 571 
Sawyer H.B., 1953, JRASC, 47, 229 

Sawyer Hogg H., 1973, Publ. David Dunlap Obs., 3, No. 6 
Simon N.R., Clement CM., 1993, ApJ, 410, 526 
Simon N.R., Teays T., 1982, ApJ, 261, 586 
Searle L., Zinn R., 1978, ApJ, 225, 357 

Smith H.A., 1995, RR Lyrae Stars. Cambridge University Press, 
Cambridge 

Stetson P.B., 1993, User's Manual for daophot ii 

Storm J., Carney B.W., Beck J.A., 1991, PASP, 103, 1264 

van Albada T.S., Baker N., 1971, ApJ, 169, 311 

VandenBerg D.A., Swcnson F.J., Rogers F.J., Iglesias C.A., 

Alexander D.R., 2000, ApJ, 532, 430 
Walker A., 1994, AJ, 108, 555 
Walker A., 1998, AJ, 116, 220 

Walker A.R., Kovacs C, 2000, in The Impact of Large-Scale Sur- 
veys on Stellar Pulsations, eds. L. Szabados & D. Kurtz (San 
Francisco: ASP), ASP Conf. Ser., 203, 268 

Welch D.L., Stetson P.B., 1993, AJ, 105, 1813 

Zinn R., West M., 1984, ApJS, 55, 45 



(http: / / www.astro.utoronto.ca/~cclemcnt/cat/clusters. htm] ) 
Clement C, Shelton I., 1999, ApJ, 515, L85 
Davis H., 1917, PASP, 29, 260 



© 2000 RAS, MNRAS 000, 000-000 



